Exchange interaction between magnetic adatoms on surfaces of noble metals 
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We present first principles calculations of the exchange interactions between magnetic impurities 
deposited on (001), (110) and (111) surfaces of Cu and Au and analyze them, in particular, in the 
asymptotic regime. For the (110) and the (111) surfaces we demonstrate that the interaction shows 
an oscillatory behavior as a function of the distance, R, of the impurities and that the amplitude of 
the oscillations decays as 1/R^. Furthermore, the frequency of the oscillations is closely related to 
the length of the Fermi vector of the surface states existing on these surfaces. Due to the asymmetry 
of the the surface states' dispersion, the frequency of the oscillations becomes also asymmetric on 
the (110) surfaces, while on the Au(lll) surface two distinct frequencies are found in the oscillations 
as a consequence of the Bychkov-Rashba splitting of the surface states. Remarkably, no long range 
oscillations of the exchange interaction are observed for the (001) surfaces where the surface states 
are unoccupied. When burying the impurities beneath the surface layer, oscillations mediated by 
the bulk states become visible. 
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I. INTRODUCTION 

The Rudermann-Kittel-Kasuya-Yoshida (RKKY) in- 
teraction in bulk materials is known for a long time,^ 
and has become a textbook knowledge. In the past two 
decades, also the Oscillatory Exchange Coupling in mag- 
netic multilayers has been extensively studied both in 
experiment^iS, and in theorj*^^— . Recently, surface nanos- 
tructures containing only a few, or even just one atom 
were fabricated, and the exchange interaction between 
the individual atoms could be measured directly^'S. Var- 
ious aspects of this interaction have already been dis- 
cussed from first principles j^^— 

In bulk, the magnitude of the RKKY interaction de- 
cays as where R is the distance between the im- 
purities, and the extremal vectors of the Fermi surface 
determine the frequency of the oscillatory interaction. It 
is well known that the (111) surface of noble metals con- 
tains Shockley type surface states^""—, which behaves as 
a two dimensional free electron gas (2DEG). It has been 
shown even earlier that the 2DEG mediates an RKKY 
type of interaction between magnetic impurities which 
has the asymptotic form, sm{2kpR)/ , where kp is the 
radius of the Fermi surface (circle) .^-"^^ 

The study of realistic, more complicated (non-free- 
electron-like) Fermi surfaces or the absence of a (par- 
tially) occupied surface state call for more elaborate the- 
oretical tools. Such methods account for the semi-infinite 
host system beneath the surface, by calculating the sur- 
face Green's function as based on first principles 
Previous studies of the exchange interaction at surfaces 
concentrated merely on the (001) and (111) surfaces of 
the host material being mostly copper jS^— and a compar- 
ative study of the vicinal surfaces of different host ma- 



terials is still missing or incomplete, especially regarding 
the asymptotic regime. 

In the present work, we perform calculations for the 
(100), (110) and (111) surfaces of Cu and Au in or- 
der to understand the role of various surface proper- 
ties and their influence on the frequency and amplitude 
of the exchange interactions between Co atoms placed 
on these surfaces. It is well-known that in case of 
Au(lll) the large spin-orbit interaction results in a split- 
ting of the surface states' dispersion, called the Bychkov- 
Rashba splittingji^iilii2ii^i2^ Here we demonstrate how 
the Bychkov-Rashba splitting manifests itself in the ex- 
change interaction on Au(lll) and also on Au(llO). 



II. METHOD OF CALCULATION 

The Screened Korringa-Kohn-Rostoker (SKKR) 
metho d^^'^" combined with the embedding technique^^ 
enables a precise treatment of a finite cluster of impu- 
rities embedded into a two-dimensional translational 
invariant semi-infinite host. Within multiple scattering 
theory, the electronic structure of a cluster of embedded 
atoms, C, is described by the corresponding scattering 
path operator (SPO) matrix, Tc(e) (e being the energy), 
given by the following Dyson equation, 



(1) 



where th{e) and Th{e) denote the single-site scattering 
matrix and the SPO matrix of the unperturbed host sites 
at the place of cluster C, respectively, while tc{e) stands 
for the single-site scattering matrix of the embedded 
atoms. Self-consistent calculations within the local den- 
sity approximation (LDA) of the density functional the- 
ory (DFT) can then be easily performedi^ In the present 
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calculations no attempt was made to include surface re- 
laxations: the cluster and the host sites refer to the posi- 
tions of an ideal FCC lattice with the experimental lattice 
constant of Cu (ofcc — 3.6147 A) and Au (afcc = 4.0648 
A). First, we performed self-consistent (relativistic) cal- 
culations for the (001), (110) and the (111) surfaces of 
Cu and Au, then for Co impurities placed on the top 
of these surfaces. The details of the scf calculations are 
described in Ref. 

Once the self-consistent potential and exchange field 
for a single impurity has been obtained, we used the mag- 
netic force theorem'^^ to calculate the interaction between 
two magnetic impurities. As in this procedure the relax- 
ation of the electronic structure due to the proximity of 
the impurities is neglected, the exchange interaction be- 
tween two impurities should be calculated from the grand 
canonical potential, f2 = E'band — £fN, with i?band, e_F 
and N being the band-energy, the Fermi level and the 
number of valence electrons, respectively. 

For uniaxial systems, such as considered in this paper, 
the dependence of the grand potential on the orientations 
of the magnetic moments, ei and 62, can be written up 
to second order as^i^^ 

(2) 

where K is the uniaxial anisotropy constant of a single- 
impurity and J(i?) is a 3 X 3 matrix comprising the 
isotropic, symmetric anisotropic and the asymmetric ex- 
change (Dzyaloshinskii-Moriya) interactions between two 
impurities placed at a relative position, R. 

As usual, we define the exchange coupling energy be- 
tween two impurities as 

J(e; R) = n{e, e; R) - n{e, -e; R) , (3) 

a quantity that still depends on the orientation, e. From 
Eq. ([2]) we, however, see that 

J(e; R) = J(i?) + e J^y™(i?) e , (4) 

where 

J{R) =TrJ{R)/3 , (5) 
is the isotropic exchange coupling parameter and 

jsym(^) = \ (J(i?) + rrJ(i?)) - J(i?)I , (6) 

is the traceless symmetric part of J(i?) with the unit 
matrix I. It should be noted that the so-called symmet- 
ric exchange (two-site anisotropy), J*'y™(^), arises purely 
from relativistic effects and it scales with the squared 
spin-orbit coupling strength. Since our calculations in- 
dicate that it is by several orders smaller in magnitude 
than the isotropic exchange coupling, J{R), in this study 
we disregard the orientational dependence of J(e; R) and 
present calculation of Eq. ^ at the specific orientation. 



e = z. In order to carefully trace the long-range oscilla- 
tions of J(i?), when calculating the host Tt matrices, see 
Eq. ([T]) , we used 30 000 fc| | points in the irreducible wedge 
of the surface Brillouin zone, while 10 000 points were 
sufficient to use in the absence of occupied surface states. 

In Table I, calculated exchange interactions are shown 
between impurities at selected distances for the case of a 
Cu(lll) surface. Note that negative/positive values indi- 
cate ferromagnetic/antiferromagnetic coupling between 
the spins. Two approaches are compared: (i) using sin- 
gle impurity potentials as described above and (ii) tak- 
ing into account the proximity of the two impurities in 
terms of self-consistent calculations. We denote the re- 
spective results for the exchange interactions by Jii{R) 
and J2i{R)- As expected, for impurities at large dis- 
tances {R — 10 a and 20 a, where a is the nearest neighbor 
distance) the two methods are in reasonable agreement. 
For the nearest neighbor pair, however, the exchange en- 
ergy obtained from self-consistent two-impurity poten- 
tials are roughly twice as big as the ones calculated from 
the single-impurity approach. 

R/a Jli Jli 

1 -197 meV -106 meV 
10 -24.5 ^eV -24.3 ^eV 
20 1.29 ^eV 1.31 ^ieV 

TABLE I: Exchange interactions between Co impurities at 
selected distances, R (in units of the nearest neighbor dis- 
tance, a), on a Cu(lll) surface as calculated by two different 
approaches (see text). 



III. RESULTS AND DISCUSSION 

A. Nearest neighbor interactions 

Although the main purpose of the current paper is to 
investigate the asymptotic behavior of the magnetic in- 
teraction between two surface impurities, we would like 
to briefly comment on the nearest neighbor interactions 
as well. For two impurities that are close enough to have 
sufficient overlap of their wave functions, it is the di- 
rect exchange which dominates the interaction. It gives a 
strong, but short range coupling which decreases rapidly. 
We compare the exchange energy for nearest neighbor 
impurities in bulk and on surface as calculated from the 
single impurity approach (see above). As we have seen, 
for the nearest neighbor pairs this is a crude approach, 
but it provides a qualitative estimate for the desired com- 
parison. 

Our results for Cu and Au hosts are presented in Ta- 
ble ini In both cases the coupling of the two Co spins 
are ferromagnetic in the bulk. This interaction is largely 
enhanced at the (100) and (111) surfaces. In general, 
this enhancement can be correlated with the decreased 
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number of host atoms in nearest neighbor position below 
the Co atoms, 4 for (100) and 3 for (111), correspond- 
ingly, with a decreased hybridization between the Co and 
host atoms. Though very small in magnitude, on top of 
the (110) surface the nearest neighbor exchange interac- 
tions become antiferromagnetic. This observation clearly 
demonstrates that the surface electronic structure can 
dramatically change the interactions between adatoms 
as compared to the bulk. 
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TABLE II: (color online) Calculated nearest neighbor ex- 
change interactions in bulk and on different surfaces of Cu 
and Au. 
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B. Asymptotic behavior 

We calculated the exchange interactions, J{R), be- 
tween two Co adatoms deposited on top of the (111), 
(110) and (100) surfaces of Cu and Au surfaces for dis- 
tances up to i? ~ 100 — 150 A, to be considered safely 
as the asymptotic region. The calculated results are 
shown in FiglT] as a function of the distance between the 
adatoms. In case of the (100) and the (111) surfaces, the 
two Co atoms were placed along the (110), i.e., the near- 
est neighbor direction, while in case of the (110) surface 
along the (001), i.e., the next nearest neighbor direction. 

For the (111) surface of Cu, see left-top panel of Fig. [U 
J{R) shows clear oscillations with a period of L = 18.5 
A, corresponding to a Fermi wavelength of the surface 
states, kp = tt/L = 0.17 1/A. This period IS very 
close to the value obtained from STM measurements and 
also from first- principles calculations by Stepanyuk et al, 
L = 15 A OT kp = 0.22 1/A."'i2 The difference between 
the theoretical results can mostly be attributed to differ- 
ent angular momentum cut-offs and to different surface 
potentials used in the calculations. In agreement with 
theoretical modelaM'—, our numerical fit also confirmed 
that the amplitude of the oscillations decays as 1/ E? . 

The asymptotic curve of J{R) displays a more compli- 
cated behavior for Au(lll), see right-top panel of Fig. [1] 
since this curve could be fitted as the sum of two oscil- 
lations with the Fermi wavelengths, k\, = 0.104 1/A and 
kp — 0.142 1/A. The appearance of the two oscillation 
periods is due to the famous Bychkov-Rashba splitting^^ 
of the Au(lll) surface states experiencing the strong 
spin-orbit interaction of Au. This splitting of the sur- 
face states gives rise to two distinct spherical Fermi-cuts, 



FIG. 1: (color online) Calculated exchange interactions as a 
function of the distance between Co adatoms on top of the 
vicinal surfaces of Cu and Au along the directions specified 
in the text. Symbols refer to the calculated data, solid lines 
to the fitted curves. 



thus, to two distinct asymptotic oscillations for the ex- 
change interaction between magnetic adatoms. It should 
be noted, however, that, mainly because of the imprecise 
treatment of the surface potential within the ASA, our 
values for the Fermi wavelengths are typically smaller, 
while their difference, Akp = 0.038 1/A, is larger as re- 
ported in the literature J^"— 

As studied by Petersen et al. in terms of Fourier 
transform STM^ there exist partially occupied Shock- 
ley states on the Cu(llO) surface that are located at 
the boundary (Y point) of the two dimensional Brillouin 
zone. These surface states naturally give rise to oscilla- 
tory exchange interactions between magnetic adatoms as 
demonstrated in the middle panels of Fig. [TJ both for 
Cu(llO) and Au(llO). 

Noticeably, the shape of the cuts of these surface states 
is elliptic due to the C2v symmetry of the surface. Con- 
sequently, the period of the oscillations depends on the 
direction along which the two adatoms are pulled apart <^ 
Fig. [2] presents the exchange interactions in case of 
Cu(llO) along the (001) and the (110) directions. The 
Fermi wavelengths obtained from fitting these curves, 
0.141 1/A and 0.172 1/A along the (001) and the (110) 
directions, agree indeed very well with the values derived 
directly from the dispersion relation of the surface states, 
0.138 1/A and 0.164 1/A respectively. 
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In case of Au(llO), we observed two oscillations in 
J{R) along the (001) direction, corresponding to the 
Fermi wavelengths, 0.098 1/A and 0.118 1/A. Similar 
to Au(lll), this can again be explained by the fact 
that the Au(llO) surface states experience an anisotropic 
Bychkov-Rashba splitting.™ Note, however, that the am- 
plitude of the oscillation with the longer period {kp — 
0.098 1/A) turned out to be about 16 times larger than 
the one for the shorter period {kp = 0.118 1/A). Most 
likely, an even larger imbalance of the amplitudes ap- 
plies to the asymptotic exchange interactions along the 
direction (110), since in this case it was not possible to 
resolve numerically the two different frequencies in the 
oscillations: our fit confirmed only the short period with 
kp = 0.151 1/A. 




50 100 150 

R{k) 

FIG. 2: (color online) Calculated exchange interaction on the 
Cu(llO) surface along two different directions: (001) (circle) 
and (110) (filled square). The solid and dashed lines corre- 
spond to the fitted curves. 



In obvious contrast to all previous cases, in case of the 
(100) surfaces (lower panels of Fig. [IJ the exchange in- 
teractions show no oscillations in the asymptotic regime. 
For Cu(lOO), J{R) is well described by an exponential 
decay. For Au(lOO), it also rapidly decreases with R, 
though, J(i?) changes sign at i? ~ 30A. Therefore, we 
may conclude that, lacking (occupied) free-electron type 
surface states, there is no long-range RKKY type ex- 
change interaction between the magnetic atoms on the 
(100) surfaces. At first glance, this statement contradicts 
to the observation of Stepanyuk et alr^ who established 
an oscillatory coupling even in case of Cu(lOO) . Their cal- 
culations were, however, restricted to distances, i? < 10 
A, that can not be regarded as the asymptotic regime. 

However, once the Co adatoms are placed in the sur- 
face layer, an oscillatory exchange interaction appears, 
although with a different frequency than in the bulk, see 
right-top panel of Fig. [3] for the case of Cu(lOO). When 
placing the impurities even deeper beneath the surface, 
the frequency of the bulk RKKY interaction is quickly re- 
covered. It is well-known that the Cu Fermi surface has 
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FIG. 3: (color online) Calculated exchange interactions be- 
tween two Co impurities placed at different vertical positions 
with respect to a Cu(lOO) surface: S -f 1 on top, S in the 
surface layer, S — 1 in the subsurface layer, as well as in the 
bulk. 



four extremal vectors along the (110) direction.- Among 
them, the largest one, related to the (110) diameter of the 
'dog-bone'-shaped Fermi surface, causes the rapid oscil- 
lations in the bulk as seen on the bottom-right panel of 
Fig. [31 Interestingly, our numerical analysis shows that 
the frequency of the oscillations for adatoms in the sur- 
face (S) layer correlates with the much shorter extremal 
vector of the Fermi-'neck'. 

So far, we presented numerical results for the exchange 
interaction between two Co impurities. In Fig. 2] the 
exchange interactions between two Co adatoms are com- 
pared with that between two Fe adatoms placed on top 
of Cu(llO) along the (111) direction. From this figure, it 
can clearly be inferred that the frequency of the asymp- 
totic oscillations is the same for both cases: the type of 
the adatoms becomes manifest in the amplitude and in 
the phase of the oscillations. The period of L = 20.4 A 
corresponds to a Fermi wavelength of kp = 0.154 1/A 
which, reassuringly, coincides well with the length of the 
spanning vector of the elliptical Fermi surface cut along 
the corresponding direction. 



IV. CONCLUSION 

We presented first principles calculations of the ex- 
change interactions between magnetic impurities de- 
posited on the vicinal surfaces of Cu and Au. In full 
agreement with previous theoretical studies, for the (110) 
and the (111) surfaces we demonstrated that the inter- 
action in the asymptotic regime is oscillatory: the am- 
plitude decays as 1/R^, where R is the distance of the 
adatoms, and the frequency of the oscillations coincides 
with the length of the Fermi vector of the Shockley type 
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FIG. 4: Fitting curve of the calculated exchange interactions 
between two Fe (dashed line) and Co (solid line) adatoms 
placed on top of a Cu(llO) surface along the (111) direction. 



This, in turn, resulted in an anisotropic periodicity of 
the exchange interaction when varying the direction be- 
tween the adatoms. Moreover, for the Au(lll) surface 
two distinct frequencies are found in the oscillations at 
any direction as a consequence of the Bychkov-Rashba 
splitting of the surface states. 

Our most remarkable observation is the lack of long 
range oscillations in the asymptotic exchange interaction 
for the (001) surfaces. This finding can be correlated 
with the fact that there are no partially occupied surface 
states in this case. We should note that the bulk Fermi- 
surface related oscillations decaying as predicted 
by Lau and Kohn^^ could not be numerically resolved 
in our calculated data. Burying, however, the impurities 
beneath the surface layer, oscillations mediated by the 
bulk states become apparent. 
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